Levy flights on the half Hne 
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We study the probability distribution function (pdf) of the position of a Levy flight of index < 
Q < 2 in presence of an absorbing wall at the origin. The solution of the associated fractional Fokker- 
Planck equation can be constructed using a perturbation scheme around the Brownian solution 
(corresponding to a = 2), as an expansion in e = 2 — a. We obtain an explicit analytical solution, 
exact at the first order in e, which allows us to conjecture the precise asymptotic behavior of this 
pdf, including the first subleading corrections, for any a. Careful numerical simulations, as well as 
an exact computation for a = 1, confirm our conjecture. 



I. INTRODUCTION 

Random walks are trajectories consisting of a collec- 
tion of random steps. They are employed to model 
the stochastic activity observed in many fields such as 
physics, biology, quantitative finance or computer sci- 
ence. As such they have been widely studied by mathe- 
maticians [iHSj and physicists [IHZ] . One of the simplest 
example is given by a symmetric one dimensional ran- 
dom walker whose position, a;(n), after n steps evolves, 
for n > 1, according to 

x{n) = x{n — 1) + 77(71) , (1) 

starting from a;(0) = 0. Moreover we consider random 
steps independent and identically distributed, according 
to a probability distribution (p{r]): consequently the ran- 
dom walk is Markovian and homogeneous. 

Despite being simple to define, most of the proper- 
ties of a random walk remain difficult to determine an- 
alytically. However, when the number of steps is large, 
the random walk displays a "universal" behavior and the 
statistics of the position x{n) becomes independent of 
most of the details of >p{ri). This asymptotic regime is 
the one we have more chance to characterize and it is 
often the one which is relevant for applications. In this 
limit, two cases should be distinguished. If ip{r]) is a 
narrow distribution the random walk exhibits a diffusive 
behavior, x{n) ~ n^/^ and only the second moment of the 
step distribution, r]'^ip{r])di] = a^, affects the statis- 
tics of x{n). On the other hand if the random variables 
77(71) 's have a broad distribution, with a diverging second 
moment, i.e. 



ll+a 



\v\ > 1 
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random walk exhibits a super- 
n^/". Such power law distribu- 



with < a < 2, the 
diffusive behavior x{n) 

tions ([2]) have been initially studied in the early sixties in 
economy [Sj and finance [9 and in the early eighties they 
started to proliferate in physics where they have found 
many applications ranging from disordered and glassy 
systems, super-diffusion in micellar systems, laser-cooling 
of cold atoms TD], random matrices [TT], photons in hot 
atomic vapours 12. etc... One striking feature of such 
processes is that their statistical behavior is governed by 
a few rare events, whose occurrence are thus governed by 
the tail of the distribution. 

Also in this case, when the number of steps becomes 
large, we expect that the statistics of x{n) becomes inde- 
pendent of the details of ip{ri) except for the index a and 
the constant c. In particular, in absence of boundaries, 
the central limit theorem ensures that the propagator, 
P{x, n), i.e. the probability to find the particle in x after 
77, steps, converges to a stable distribution given by: 



P{x, n) 



Riy) , 



^''y rill — p^\°-^\° 



(3) 
(4) 



where the parameter a is related to the constant c [see 
Eq. ([s]) below] and where we have introduced, from di- 
mensional analysis, the rescaled variable 



y ■ 



(5) 



Although the Fourier transform of R{y), its characteristic 
function, has a very simple expression (|4]), there is no 
simple closed form expression for R{y), except for a = I 
which corresponds to Cauchy distribution. It however 
admits the following asymptotic expansion, valid for any 
value of a [3: 



Riy) 



1 

7r|y| 



k=l 



k+l 



kl\y 



ak 



• (6) 



2 



One thus sees on this expression that R{y) inherits the 
power law tail of the step distribution ^ . One can fur- 
ther show [U |3] that the amplitude itself is not renormal- 
ized such that to leading order 



R{y) 



,l+2a 



0{y 



-1-3q 



which fixes the value of i 



3in(^f )r(a + l) ' 



(7) 



(8) 



while d can also be obtained explicitly from Eq. (|6|. 

Although free Levy flights are thus perfectly well un- 
derstood, there are physical situations which actually in- 
volve Levy flights in a confined geometry. An interesting 
example is the Levy flight model which has been proposed 
[13] to describe the transports of solar photons in cloudy 
atmosphere. These photons are eventually reflected back 
to space or absorbed by the ground, so their trajectories 
are bounded random walks. In such cloudy atmosphere, 
whose heigth is typically of the order of 10 kms, the pho- 
tons can be trapped in optically dense region (inside the 
clouds), travelling less than a meter between scatterings, 
while they can "fly" many kilometers from cloud to cloud. 
It was shown experimentally that Levy flights provide a 
reliable description of the photons transport in such sit- 
uations |14j . It was shown, in addition, that the one- 
dimensional model is a reasonable approximation of the 
three-dimensional geometry |15) . More recently. Levy 
flights in conflned geometry have also found applications 
in the context of random search problems |16j . 

Obviously, when the walker is confined inside a do- 
main, the central limit theorem does not apply, but the 
scaling analysis is still valid and the "universal" behav- 
ior of the rescaled position y = xjn}-!"- is still expected. 
Computing the statistics of the rescaled variable y in 
presence of confinement is in general possible for Brown- 
ian motion (a = 2) for which powerful analytical tools are 
available such as path- integral techniques [171 HH] • Un- 
fortunately, for Levy flights, analytical approaches are 
usually quite difficult. Recently a lot of papers adver- 
tised the possibility to write a Fractional Fokker-Planck 
equation for Levy flight propagator [6] : 

^P(x,t) = a"^P(x,t) , - 0) = <5(x) , (9) 

where the continuum time t captures the large n behav- 
ior of the random walk and the fractional operator, gj^, 
is the Riesz-Feller derivative of fractional order a > 
[19|, I20j , which has an integral representation involving a 
singular kernel of power-law form. In absence of bound- 
aries this equation can be simply written in Fourier space 



di 



P{k,t) = -\ak\°'P{k,t) , P{k,t = 0) = 1 , (10) 



and it is easy to check that the free propagator introduced 
m Eq. (H is a solution of Eq. Um, (with the identifica- 
tion, t ^ 1). 



In presence of boundaries the translational invariance 
is broken and the Fourier representation becomes use- 
less. For a — 2, when the fractional operator becomes 
the standard Laplacian, the method of images allows to 
express the propagator in presence of boundaries as a 
linear combination of free propagators. Unfortunately 
we will see that these techniques can not be applied, for 
a < 2 [51] [22] . More generally, if translational invariance 
is lost, the fractional Fokker-Planck equation in Eq. (|9|, 
becomes a difficult integro-differential equation with non- 
local boundary conditions. One could conclude that the 
Fokker-Planck formalism is of little help for Levy fiights. 
However, following a recent work by Zoia, Rosso and Kar- 
dar [23] , we show that Eq. (|9] ) can be studied also in 
presence of boundaries using a perturbation theory where 
the small parameter e <C 1 is a = 2 — e. At variance with 
the method of Ref. [23] we perform this perturbation the- 
ory directly in the continuum limit without resorting to a 
discretization (in space) of the trajectories. The calcula- 
tions turn out to be somewhat simpler in this continuum 
setting. 



Final position- 



FIG. 1. Trajectory of a random walker with an absorbing 
boundary condition at the origin. P+{x,t), computed pertur- 
batively in Eq. ( |C18[ ), is the probability density function of 
the position of such a random walker at time t. 



For concreteness we will study in detail the case where 
there is an absorbing wall in a; = as depicted in Fig. [T] 
we thus consider only the paths that remain positive up 
to the n"^ step. In the limit of large n the probability 
density function to find the particle in x after n time 
steps also takes the scaling form, as in Eq. (|3]): 



P+ix,n) 



1 



7,1/" 



R+ iy) 



(11) 



We emphasize that R+{y), being a probability density 
function, is normalized, i.e. R-^(y)dy — 1, while, in 
Ref. [21], Zumofcn and Klafter studied a similar quan- 
tity, which is however not normalized (see also Ref. 24 ) . 
In particular, they were able to show that the small ar- 
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gument behavior of i?+ (y) is given by [STJ [5S] 

R+{y) - 2/"/' , (12) 

in contrast with the method of images that would pre- 
dict i?+(y) ^ y. Our perturbative approach allows us to 
conjecture the exact behavior of the tail of R+{y), which 
controls the statistics of rare events : 



R+iv) = 4t- 



o(y-2-"), 2>a>l, 

(13) 

\-o{y-^-^°') , 1 > a > . 



To our knowledge, only the exponent of the leading term 
R+{y) oc y~^~", was known from Ref. [H]. Here we 
obtain the exact result for the amplitude c+ 



c+ = 2c 



(14) 



where c is the amplitude of the tail of the jump distri- 
bution ( (2I). In addition, the first sub-leading corrections 
in Eq. (|13|), by comparison to the free case ([7|, also 
bear the fingerprints of the absorbing wall. These results 
(131 are first obtained analytically for a close to 2, i.e. 
a = 2 — e, using a perturbation theory to first order in e. 
We also obtain this behavior for a = 1, albeit with log- 
arithmic corrections for the subleading term, for which 
an exact calculation can be done. We then demonstrate 
this behavior using thorough numerical simulations. 

The paper is organized as follows. In Section II 
we present the general framework of the perturbation 
scheme. We first illustrate it on the simplest example of 
the free propagator in Section II A, and then in Section 
II B we study the case with an absorbing boundary at 
the origin. The discussion of the results is left in Section 
II C. Section III contains the results of our numerical 
simulations and our conclusions are in section IV. The 
Appendices A, B, C contain some technical details. 



II. PERTURBATION SCHEME 

In this section, we set a — 1 for simplicity and without 
loss of generality, and write the Fractional Fokker-Planck 
equation (|9| in the familiar Schrodinger form 



dtP{x, xq, t) = UPix, xo, t) , 
P{x, Xq, i = 0) = 5{x — Xo) , 



(15) 
(16) 



where the propagator P(a;, xq, t = Q) represents the prob- 
ability density to find a particle in the interval [x, x -\- dx\ 
at time t, knowing that the particle was in xq at time 
and the operator V. is the fractional operator of index 
a. In Quantum Mechanics Eq. ( 15 ) corresponds to the 



Schrodinger equation of the element {x,xq) of the density 
matrix P at the temperature 1/t. The general solution 
of Eq. ( 15 1 reads 



P(x,a;o,t)= / d(?V:(2:o)V'9(a:)e^(«^* 



(17) 



where E{q) are the eigenvalues and ipq{x) the associated 
eigenfunctions of the operator H. They are solutions 
of the eigenvalue problem, Hipqix) — E{q)ipq{x), with 
the appropriate boundary conditions and satisfy the or- 
thonormality and closure relations: 

J^q{x)rq'{x)dx = S{q-q'), (18) 

4>qix)^p*g {x')dq = S{x - x') . (19) 



Please note that the domain of integration over x depend 
on the boundary condition and the integration over q is 
meant over the all spectrum. If the spectrum is discrete 
the integral is replaced by a discrete sum. 

The eigenvalue problem of the operator H is in general 
very difficult and the solution is known only in the ab- 
sence of boundaries where the eigenfunctions are simple 
plane waves. The case a — 2 corresponds to the standard 
Laplacian for which a lot of results are known. In partic- 
ular the spectrum can be exactly solved in one dimension. 
A possible strategy to gain insights on the behavior of a 
Levy flights in presence of boundaries is to use the per- 
turbation scheme around the well studied a = 2 case. 
For a = 2 — e the operator H and the propagator can be 
expanded in powers of e, 

P{x,xo,t) = P^"\x,xo,t) ~ eP^'\x,xo,t) + Oie^) 

n^no-eni + o(e^) , (20) 

where P^^\x, xq, t) is the propagator associated to Hq = 
with the prescribed boundary conditions. T he e x- 



pression of the first correction P^^' (x, xq , i ) in Eq. ( 20 ) is 



well known from Quantum Mechanics and it is re-derived 
in Appendix [a| It reads, writing P^^^ = P*^^-' (x, xq, i): 



P 



(1) _ 



^E{q)t _ ^E(q')t 
^^JqixW.ix,) ^(^,)_^(^) 



q\niW), 



. . (21) 

where we use the notation = J dq and where it is un- 
derstood here that ipq{x) are the eigenvectors and E{q) 
the corresponding eigenvalues of Hq and for the matrix 
elements {qlHilq') we use the bra-ket notation, borrowed 
from Quantum Mechanics with {x\q) = iliq{x). This for- 



mula ( 21 ) is the cornerstone of the perturbation approach 



presented here. 



Absence of boundaries 



We first illustrate this perturbative approach by com- 
puting the first order correction to the Brownian propa- 
gator in absence of boundaries, i. e. the Gaussian prop- 
agator, when a = 2 — e and xq — 0. In this simple case 
the eigenvalue problem of Ho = &i gives: 



1 



'2tt 



e^'\ E{q) 



-oo < q < oo. (22) 



4 



Using that fc^ = fc^ — e fc^ log |fc| + O(e^), the matrix 
element {q\'Hi\q') can be explicitly computed: 

{q\'^Li\q') = JJJ dxidx2'^ip*{xi)i)q,{x2) 

X fc2 log |fc|e*'=(^i"^^) = (5(q - q')q^ log |q|, (23) 

where the integral over xi,X2,k is performed over the 
whole real axis. Note that, thanks to the oscillating term 
^ik{xi-x2) ^ the integral over k in (23) is dominated by 
the small values of k where one can safely expand fc^^*^ 
in powers of e. From Eq. (21 ) one obtains: 



dq 

P''^\x,t) — —t / — cosgxq^ logge" 
Jo ^ 



Using k = q^/i and the scaling variable z 
Eq. ( 24 ) can be recast in a simpler form 



yrtP^^^\x,t) ^ Ra{z) + RB{z)\ogt , 

f°° dk 2 
Ra{z) — — I — cos kzk^ \ogke~^ 

Jo 

Rb{z) ^ ( cos kzk^ e'^"" 

Jq 



27r 



From scaling argument one expects ^ 



P{x,t) 



-R 



(24) 

x/y/i, 

(25) 
(26) 

(27) 
(28) 



and a general issue of these perturbative computations is 
that for a = 2 the natural scaling variable is z = x/ y/i, 
while for a = 2 — e the correct scaling variable is y = 
x/t^^ , which also admits a perturbative expansion. For 
this reason our final result should be recast in terms of 
in order to identify the perturbative expansion of the 
scahng function R{y) in Eq. ( 28 ) as 



i?(y) = i?W(y)-ei?(i)(y)+0(e2), (29) 
where i?^°^(?/) is the Gaussian propagator given by 

1 



2y^ 



(30) 



This can be done if in the equation t^^°'P{x, t) — R{y) we 
expand at the first order in e both i^/" ~ -y/i + |-\/ilogt 



and y 
write 



\ez\ogt. After some simple algebra we can 



^rtP'^^\x,t) = R^^Hz) + -\ogt{l + zd,)R'^'^\z) . (31) 



Comparing with Eq. (25) we identify 4i?B(z) = i?(°)(z) + 
zd^R'^^\z) so that we conclude that R'^'^>{y) in Eq. ^ 
is given by: 



R'-^\y) = - / — coskyk^ log/ce" 



(32) 



By performing an asymptotic analysis of Eq. ( 32 1 for 



Eq. (29) which converges nowhere but exists as a for- 



mal power series: 

(2fc)! 1 
•^^"^{k - l)\ y'^''+^ ^ ^ \y^ ' y 



1 12 180 



5 y7 



This result in in perfect agreement with the expansion 
given in Eq. ^ for a = 2 — e. 



B. Propagator in presence of an absorbing 
boundary at a; = 

We consider now G-f.{x,xo,t) dx, the probability to 
find a particle in the interval [x, x + dx] at time t, know- 
ing that the particle was in xq at time and given that it 
stayed positive up to time t (see Fig. [ij. At variance with 
the free propagator, the integral over x of G+(x, xo,t) 
is smaller than one and gives the fraction of surviving 
walker up to time t {i.e. the survival probability). In the 
geometry defined by Eq. ([T]) the initial position of the 
discrete random walk is xq = 0. Here we are consider- 
ing a process which is continuous in time and this initial 
condition xq = 0, together with the presence of an ab- 
sorbing boundary at the origin is ill-defined. Indeed, it is 
well known that if the continuous time walker crosses zero 
once it will re-cross zero infinitely many times immedi- 
ately after the first crossing. Therefore, it is impossible to 
enforce the constraint a;o = and simultaneously forcing 
the position of the continuous time walker to be strictly 
positive immediately after. Therefore we set xq > and 
small in order to regularize the continuous time process 
and we will take the limit — >• at the end of the cal- 
culation. Our final result corresponds to the geometry of 
the Eq. ([T]) in the limit of a large number of steps. 

It is useful to express G+{x,XQ,t) in terms of rescaled 
variables as in Eq. ([5]): 



G+{x,xo,t) = ■Jj^Z{y,yo), 



(33) 



The scaling function Z{y^ yo) depends explicitly on a and 
we compute it here in perturbation theory for e = 2 — 
a < 1: 



Z{y,yo) = Z^"\y, 



Vq) 



^-Z'^'Hy, 



Vo) 



0{e^). (34) 



In presence of an absorbing boundary at the origin, the 
action of the fractional operator can be written as 

j dx'U^')j^ _(_|fc|«)e'M--')=i?(g)V;,(a;), 

(35) 

and the solution is known only for a = 2: 



large y one finds a series expansion of R{y) given in 



ijq{x) = 0{x)^^ Sin(gx) , E{q) = -q' , q > , (36) 

where 9{x) is the Heaviside function: d{x) = 1 if 
X > and 9{x) = if cc < 0. At zeroth order in e. 
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the scaling function Z^^^{y,yo) can be computed from 
Eq. (171 with ipq{x) given in Eq. (36). Using the identity 
2 sin(^) sm{kyo) = cos [k{y - yo)\~ cos [k{y + j/q)] one 
obtains: 



Z^'^\y,yo)^Oiy)0iyo) 



dk 



e ''^ {cos[k{y - yo)] 



- cos[k{y + yo)] 



(37) 



Let us note that the same result can be straightforwardly 
obtained using the method of images: 

Z^'\y,yo) = e{y)9{y,)[R^°Hy-yo)-R('^Hy+y,)] . (38) 

At the first order in e we first compute the matrix el- 



ement {qlHilq') which has the form given in Eq. (23) 



with the prescription that the integrals over xi and X2 
are performed over the interval (0,cxi) and tpqi^) sre the 
eigenvectors given in Eq. (36 1. The integrals over xi and 



X2 need to be regularized to be well defined and this can 
be done via the identity: 



lim 



dx 



sm{qx) — PV 



q 



TT{k^ 



+ -{6{q-k)-6{q + k)) , (39) 

where PV indicates a principal value. After some algebra, 
left in the Appendix [b| one gets [see Eq. (Bll)]: 



{q\ni\q')=S{q-q')q'log\q\ + 



2(g + <z')' 



(40) 



Combining the latter equation with Eq. ( 21 ) we can write 

an expression for p'^\x,XQ,t). Analogously to the case 
of the propagator in absence of boundaries, the inte- 
grals involved in the expression of p'^\x,Xf),t) can be 
naturally recast in term of the variables z = x/ \fi and 
zo = xqI \/t instead of the correct scaling variables y and 
yo- Following the same lines of the previous discussion 
we easily write the scaling function Z*-^-* (y, i/o) as: 

(y, yQ) = ZA [y, yo) + ZB {y, yo) , (41) 
ZA{y, yo) = R'-'^ (y - yo) - R^'^ (yo + y) , (42) 



Zsiy.yo) = \ I dki dk 



■ TT{ki + k2){k^ - kf) 



X sin(fci2/) sin(fc2yo) , 



(43) 



where the integrals over ki , fc2 run over the interval 
(0,00). It is worth to stress that the term Z^ corre- 
sponds to the images method prediction while the term 
Zb represents the violation of the images prediction at 
the first order level. 

It is easy to realize that the probability density func- 
tion R+{y) is simply related to Z(y,yo) in the following 
way 



R+iy) = lim 



Z{y,yo) 



yo^oj^ dy'Z{y',yo) 



(44) 



For a = 2 one has from Eq. ( |38[ ) in the limit yo — >■ 0: 

Zi^\y,yo) = yoZ('\y)+0{yl) , Z^'Hy) = -^e- 



which yields 



R 



(0) 



(?/)=2^" 



(45) 



(46) 



The integrals in Eqs. (|42[) and ( 43 ) which give the term 



R^^\y) have to be discussed carefully and the details are 
left in the Appendix [Cj The net result of this analysis is 
that R+{y) can be written as 

R+{y) = Rf{y) (l + eW+{y) + 0{e')) , (47) 

where W^{y) can be expressed in terms of elementary 
and special functions, see Eq. ( C18| ). From this expres- 
sion, one obtains the asymptotic behaviors of i?+ (y) 
the limit y — > 0, one finds 

y 

2 4 

where k = 2 — log 2 — I7B. In particular, the small y 
behavior in Eq. (48) is consistent with i?+(y) ^ y^ in 
agreement with previous findings [211 123] . For y — >■ 00 
one finds 

2 30F ^ 32 



R+{y) 



-(logy 



0(y'logy,e' 



In 



(48) 



R+iy) ^ e 



(49) 



, y 

where the leading term, vanishing as 1/y^, is expected 
from previous analysis [21]. We notice that the right tail 
of the R+{y) has the same behavior as the right tail of 
R{y). Quite interestingly, our perturbative result shows 
that 

^ _ R+iy) 

c 



lim 



Riy) 



2 + 0(e) 



(50) 



where c and c+ are defined in Eq. (0 and Eq. ^ 
respectively. Another signature of the boundary, re- 
vealed by this perturbative calculation, appears in the 
sub-leadi ng c orrection which vanishes as 1/y^ for R+{y ) 
[see Eq. (|49|] instead of 1/y^, as for R{y) [see Eq. ((33)]. 



C. Discussion and conjectures 

It is interesting to compare our perturbative results, 
valid in principle for 2 — a <Si 1, with the exact re- 
sults which we can obtain for the special case a = 1. 
In this case, corresponding to Cauchy random variables, 
i.e. (fiirj) = Tr~^{l + rf')~^ , one can use the results which 
were obtained by Darling [26] and Nevzorov [27] in the 
context of the extreme statistics of such Levy statistics, 
to obtain an exact result for R+iz) in terms of a single 
integral: 



R+{z)^-^z j g(^^)v~^l\l~v)-^'^dv 



li^) = 



dz 



-V 
1 



(51) 



7r(l + z2)3/4 



exp 



logM 

TT 7q 1+1*2 



du 
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Its asymptotic behaviors are given by 



R+iz) 
R+{z) 



1 

2 

2 1 



2 , Z -> 

16 logz 



(52) 
oo . (53) 



On the other hand, one has in this case Q i?(y) ~ 
l/(7ry^) when y oo such that one obtains also in this 
case a ^ 1 



c 



= hm 



i?+(y) 
R{y) 



= 2 



(54) 



Based on the perturbative result obtained above ( 50 1 and 



on this exact result (|54|) we conjecture that this relation 

C4 



2c actually holds for all value of a, as stated in the 
introduction ( 14 ) . This conjecture is corroborated below 



by our numerical simulations. 

Besides, we interpret the exponent of the sub-leading 
correction in i?+(y), which decays as in Eq. (49), 
as 4 = 2 + a, with a — 2 + 0{e) while the sub-leading 
corrections in R{y) decay as 1/y^ ^ with, instead 5 — 
1 + 2a, for a = 2-1-0(6). This leads us to conjecture that 
the subleading corrections behave actually differently for 
a > 1 or a < 1, as announced in the introduction (13 1 



R+{y) 



o(y-2-"), 2>a>l, 

(55) 

|-o(?/-i-2") , 1 > a > . 



Hence a = 1 appears as a critical value regarding these 
sub-leading corrections, for which it is not surprising to 
observe logarithmic corrections ( 53 ) . This also implies 
that the coefhcient above (551 is diverging when a — > 



1 from above. Below, we will test this behavior by means 
of numerical simulations. 
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FIG. 2. Case a = 1.5. Behavior of the right tail for (/^(ry) 
(Black) and the the rescaled final position, y, of free random 
walks of n = 1000 step (Red). Histograms are performed 
using lO'^ samples. 
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III. NUMERICAL SIMULATIONS 

We consider the case where the increments distribution 
4>{7]) is the symmetric Pareto distribution: 



0(^) = 



,for |?7| > 2< 



H 

, otherwise 



(56) 



This distribution can be sampled efhciently using random 
number drawn from a uniform distribution: 



[ran(0,i)] ° 
- [ran(0,i)]" 



, with probability 



, with probability 



(57) 



where ran(0, |) is a random number in the interval (0, \). 
We construct a large number of random walks, for each 
random walk we record the final position x{n) after n 
steps and compute the correspondent rescaled variable 
y = x„/n^/". We first present our data for a — 1.5, 
for which 0(ry) — l.b/\rif'/'^ for [t^I > 2^/^. For large 
n the distribution of y should converge to the a stable 



FIG. 3. Case a = 1.5. Free random walks of n = 1000 steps 
(Red) and random walks constrained to be positive (Blue). 
Histogram are performed using 10* samples. 



distribution centered around and with an asymptotic 
tail 



R{y) 



3 24 , _4, 



2y5/2 ' yi 



(58) 



This prediction is confirmed by our direct simulation, in 
Fig. [2] we show that the tail of 4>{'q) and R(y) coincide 
when rj^y 00. The symmetric distribution R{y) is 
also plotted in Fig. |3] where we also show (y) , the 
histogram of the rescaled final position of the random 
walks constrained to be positive. R+{y) is clearly defined 
only for positive y, vanishes at y = and, when y — >■ 00, 
decays as c+/y^/^. One of our main prediction is that, 
for large y, R+{y) / R{y) — 2 which means c_|_ = 3 for our 
model. This is verified in Fig|4] 
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FIG. 4. Case a = 1.5. Free random walks of n = 1000 steps 
(Red) and random walks constrained to be positive (Blue). 
Study of the tails. Histogram are performed using 10** sam- 
ples. The expected tails: l.^/y^^'^ and are also drawn 
(Solid Line). 




FIG. 6. Case a = 1.5. Finite size effects for random walks 
constrained to be positive. Long random walks (n — 5000 
-Blue) vs short ranomd walks (n = 250 - Red). The constant 
c+ = 3 is recovered. 




0.1 1 10 

y 



FIG. 5. Case a = 1.5. Finite size effects for random walks 
constrained to be positive. Long random walks (n — 5000 
-Blue) vs short ranomd walks (n = 250 - Red). Histograms 
are perfomed over lO'^ samples. The slope is plotted as a 
guide of the eye (Black line). 



A. Finite size effect and different values of a 

Stable distributions and universal behavior are ex- 
pected in the limit of a large number of steps {i.e. 
n — )■ oo). In our numerical simulation the asymptotic 
behavior of R{y) and R+{y) is studied for n = 1000. Is 
this number enough to capture universality? In Fig. [5] 
and Fig. [6] we study how the finite number of steps af- 
fect the function R+{y). Finite size effects are visible 
close to the boundary y — where, only for very large 
size, the distribution vanishes with the predicted expo- 
nent a/2 = 0.75. For y ~ 10 the convergence with the 



size n becomes faster and the constant 

c+ = hm R+iy)y"+^ (59) 

1/-5-00 

can be correctly estimated even with a moderate number 
of steps (see Fig|6|. 

Finally we have checked that our result for c+ apply 
to all range of < a < 2 for symmetric Levy flights. 
The asymptotic tail is more and more pronounced as 
a <^ 2. This means that the insights given by our pertur- 
bative calculation are actually valid for all Levy flights 
(see FigjT]). 




y 

FIG. 7. Random walks constrained to be positive for 
n = 1000 steps and a = 1.75,1.25,0.75,0.5,0.25 (from left 
to right). Histograms are perfomed over 10^ samples. The 
asymptotic behavior c+/y°'^^ is also plotted (Black line). 
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B. Subleading corrections 

We also check the behavior of the sub-leading correc- 
tion ( 55 ) . For a > 1 this correction is expected to behave 
like d^y'^"^^, while for a < 1 we expect that it decays 
as ^ c?+2/^^"^^. This prediction is confirmed in Fig. [s] 



for a > 1 and Fig. |9] for a < 1. 
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FIG. 8. Sub-leading correction for a = 1.25 > 1. Random 
walks constrained to be positive for n = 1000 steps. His- 
tograms are performed over 10* samples. The plateau reached 
at large y gives a numerical estimation of d+. 
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FIG. 9. Sub-leading correction for a — 0.5 < 1. Random 
walks constrained to be positive for n = 1000 steps. His- 
tograms are performed over lO'^ samples. The plateau reached 
at large y gives a numerical estimation of d+. 



This approach, following the work of Zoia, Rosso and 
Kardar [23], amounts to construct a perturbative solu- 
tion of the fractional Fokker-Planck equation (|9| with 
appropriate (non-local) boundary conditions. Here, at 
variance with Rcf. [23, the perturbation theory is car- 
ried out directly for a process which is continuous both 
in space and time. 

We have then used this perturbative method to com- 
pute, to order 0(e), the probability density function, 
of the position of such a walker with an absorbing wall 
at the origin. A different perturbation scheme (based 
on path integral) was used recently to compute the 
same quantity for the fractional Brownian motion, a non- 
Markovian process displaying anomalous diffusion. Our 
main result here is to give a precise conjecture, valid for 
any value of a, on the relation between the tail of this 
distribution and the tail of the steps of the random walk. 
Numerical simulation confirm our conjecture. 

This perturbative scheme opens the way to an analyt- 
ical study and can be used for any confined domain for 
which the Brownian solution, corresponding to a = 2, 
is known. More realistic confined geometry, like the one 
relevant to the scattering of solar photons, can be studied 
along these lines. Here we have proposed the simplest of 
these geometries. A first extension of the present study 
concerns the study of the extreme statistics of a Levy 
bridge, which is a Levy random walk on the time in- 
terval [0, t] constrained to start and end at the origin. 
Such a constrained Levy random walk, for which little is 
known, has recently received some attention in statistical 
physics |29) - in relation with some real-space condensa- 
tion phenomena -, as well as in finance |30) . Another 
interesting application of such a perturbative calculation 
could be the study of non-intersecting Levy walkers, the 
so-called "vicious" Levy walkers, which were recently in- 
troduced in Ref. [31j . 



ACKNOWLEDGMENTS 

This work was supported by the France-Argentina 
MINCYT-ECOS A08E03. R. G.-G. acknowledges the 
hospitahty at LPT and LPTMS in Orsay, A. R. and G. 
S acknowledge the hospitality at the Centro Atomico in 
Bariloche where part of this work was done. A. R. ac- 
knowledges support by ANR grant 09-BLAN-0097-02. G. 
S. acknowledges support by ANR grant 2011-BS04-013- 
01 WALKMAT. 



IV. CONCLUSION 



Appendix A: Perturbative scheme to order 0{e) 



To conclude, we have presented a perturbative ap- 
proach to the study of a Levy flight, of index < a < 2 on 
a half-line, where the perturbative parameter is e = 2 — a. 



In this appendix we give a short derivatio n o f the ex- 
pression for the first correction P*-^' (x, XQ,t) ( 20 1 given in 
Eq. (21 1. The equation for P^-^\x, xo,t) reads, see Eqs. 
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(15 20): 



where 



(Al) 

with the initial conditions P'-^'> (x, Xn .t = 0) = and 
P'-°^x,xo ,t) is given in ( |l7[ ). Eq. ( |Al[ ) is inhomoge- 
neous and we can obtain its solution using the method of 
variation of constants. We look for a solution under the 
form 

P^'Hx,Xo,t) = / c«(t)V',(x)e^(«)*dg , (A2) 



where "0? and E{q) are the eigenvectors and eigenvalues 

of 7^0 7 while the coeffi cient Cq^^(t) depend s ex plicitly on 
time t. Plugging Eqs. (A2| and (17) into (Al) one finds 



dq'e^(«')V9'(a;)5tcW(i) (A3) 
= / dg'e^(«')*V;,(xo)Hi^/','(a:) . (A4) 



We multiply the left and right hand sides of Eq. ( A3 ) 
by il'qi^) ctiid we then integrate both sides over x. Using 
the normalization relation J dxtjj*{x)ijjq'{x) — 6{q — q') 
we obtain 

e^(')*9tcW(t) = / d9'(g|Hi|g>:,(xo)e^(«')* , (A5) 



where we have used the bra-ket notation, borrowed from 
Quantum Mechanics, and such that {x\q) = ipq(x). This 
can be straightforwardly integrated, using the boundary 
condition at time t = 



cW(t)= J dq'{q\n^\q')rA^o rE^q^_E{q') 



(^E{q')t _ ^E(q)t 



(A6) 



Combining ( A6 ) and ( A2 ) one obtains Eq. ( 20 ) given in 
the text. 



Appendix B: Evaluation of an integral 

Here we compute the principal value of the integral 
entering the matrix element in Eq. ( 40 1 : 



/oo 
dk- 
-oo ^ 



qq'k"^ ln|fc| 



_ 2qq' f°° fc^lnfc 

" ^ X {k^-q^W-q'') ■ ^ ^ 

Without any loss of generality, since I{q,q') — I{q',q), 
we assume q' > q. Performing the change of variable 
k = sq, one obtains straightforwardly 



poo 

A(.)=pvy^ ds^^ 

Mx) ^ FY fds^^ 



)(s2-x2) ' 



s In s 



)(,S2 



We notice that f2{x) can be written as 



(B3) 
(B4) 



ln(s/a;) 



/.(.) = A(.)lnx + PVy^ dss^^^^^^j^^y 

(B5) 

The purpose of this trick (B5| is that the integrand has 



now only a simple pole in s = 1 so that, taking the prin- 
cipal value, yields a perfectly smooth function of x. This 
integral can then be evaluated to yield: 



f2ix) = fi{x)lnx + 



4(1 + x) ' 



so that one has 



2(7' fq'\ , , qq' 

I[q,q) = —fi - In? + TTt — TTY 
TT^ \qj 2{q + q') 



(B6) 



(B7) 



Looking at the expression of fi{x) above (B3), one sees 
that for X = 1, the integrand has a double pole in s = 1 
and one thus expects fi{x) to be highly singular in a; = 1. 
To characterize it, we compute its Fourier transform 



flip) 



dxfi{x) cos (px) 



(B8) 



which yields after straightforward manipulations of Eq. 
(B3) (which can be done e.g. with Mathcmatica) 



hip) 



cosp . 



which yields 



.1 



and finally 
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1 



5{q^q') 



Iiq,q')^Siq-q')-q^lnq+^^ 



(B9) 



(BIO) 



(Bll) 



which is used in the Eq. ( 40 1 in the text 



Appendix C: Details about the perturbative 
calculation of the propagator with an absorbing wall 

This appendix is devoted to the analysis of Z'^^^ {y, i/q) 
given by the sum of the two terms in Eqs. (41 42 43) 



in the limit yo 0. The first term Z^(y,yo)is easy to 
analyse and yields 

ZAiy,yo) = -2yodyR'^'\y) + ©(yg) , (CI) 

where i?'-^-'(j/) admits the integral representation given in 
Eq. ([32l). 
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The analysis of the small j/o behavior of Zsiy, yo) given 
in Eq. (43 1 is more subtle. To deal with this double inte- 



gral over ki and k2 , we first make the change of variable 
ki — u' and fc2 = uu' and observe that the integral over 
u' can then be performed to yield 



The integral in Xii(y,yo) contains a logarithmic sin- 
gularity when 2/0 ~^ 1 which is a bit tricky to ex- 
tract. To do so, we first perform a change of vari- 
able s = you and then add and substract the term 
yo dsie-y^/^ - e-y^/'^-y')/{s + yof. Now using 



ZB{y,yo) 



1 



du 



4y/^Jo (m-1)(u + 1)2 



ds 



1 



(s + yo) 



(e - l) = t/(log yo + log y + 7b) 



— e H I e 4,i2 _ e 



O(yolog2/o), (CIO) 



(C2) 

one obtains, for j/q — >■ 0: 



Using now the identity 



2^11 (y, yo) = (yo log yo) ye "-i + yoye ^ (7B log y) 



1 



1 



yoe 



Q(y)+0(yg log yo) 



(Cll) 



(u-l)(w + l)2 (m + 1)2 (m-1)(u + 1) 

we split ZB{y,yo) into two parts 
1 



2 XC3) 



where 

Q(y) 



^ fe^?'«-l-|-2e-Tsinh^^'' 



Zsiy^yo) 



(Ii(y,yo) +l2(y,yo)) , (C4) Finally combining Eqs. (C7), ([C9| and (Cll ) one obtains 



(C12) 



where 

2^i(y,yo) 



the small yo behavior of Zsiy, yo) in Eq. ( |C4| as 

ye" 



(«+!)= 



_ (y+^yp) _ (y-^yp) 

e " — e * 



1 / _(y-"yo)^ (y + ^yp)^ 



ZB{y,yo) 
fQ{y) 

yo 



. , (yolnyo) 
4v^ 

IE - In (y) - A{y) 



(C13) 



(C5) 



l2(^, 



(u-l)(« + l)2 



(y + ^yp) _ (y-^yp)^ 

e 4 — e " 



1 



\ I _ (y-yyp)^ 

e 



("+»'yp)" 



(C6) 



In l2(y-yo) (C6| the small yo limit can be taken easily. 
It yields 



V y 

where A(y) is given by 

(C14) 

One can then obtain Z{y,yo) = Z^^'{y,yo) — 
e{ZAiy,yo) + ZB(y, yo)) i n the limit yo ^ from the 
Eqs. (pSl), (tcTl) and ((ClSl) as: 



C'(yo Inyo) 



.-^(^-i)_ 



i2(y,yo) = yyo 



(U- 1)(U+ 1)2 \U 



1 



^(y,yo) = yo(l - ^lnyo);7^e 

Z zy TT 



ej/oZi(y), (C15) 



+ 0{yl) ■ (C7) 
We now decompose Ii{y, yo) \C6\ into two parts 

(C8) 



where Zi{z) can be read off straightfowardly from 
Eqs. (fcTl) and (ICTsl): 



2^i(y,yo) =2:ii(y,yo) +^i2{y,yo) 

, f du 

Mi{y,yQ 



Z,{y) = -2dyR^'\y) 

-T^^ 7£;-lny-^y) 

4v7r V y 



(C16) 



2^12 (y,yo) 



(^ + 1)^ 

du 



e 4 — e 



"("+1)2 



(y-yyp)^ 



(y + yyp)^ 



In Ii2(y,yo), it is straightforward to obtain the small yo 
behavior as 

ii2{y,yo) = yyo ^r^rr^e-i^i^ + O(yo') . (C9) 

Jo U [U -\- L) 



This perturbative expansion (C15) is fully consistent 
with the expected behavior ^2lt f25] 

Z{y, yo) - y! Z{y) , yo ^ , Z{y) - {y) - eZ^) (y) , 

(C17) 



where Z^^\y) is given in Eq. (C16). These integrals that 
enter the definition of Z^^\y) can then be computed, 
for instance using Mathematica, to yield the following 
expression of R+{y) from Eq. (|44]): 
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R+{y) = Rfiy){l + eW+{y) + 0{e')) , 
1 



(C18) 



60 - 2AjE - 48 log 2 - y^(18 - 67^ - 12 log 2) - 6(2 - y^y)e' 



W+{y) 



+y'{y -4)^ 2F2{ 1,1;-, 3 



24^71 



-Erf (I) - 37r(y2 _ 2)Erfi (|) + 3{y^ ~ 4)Ei 



where i?^'' (y) is given in Eq. 
a hypergeometric series |:32 . 



(|46|) and 2F2 {1,1; l,3;u) is 



From this expression (C18) 



it is straightforward to obtain the asymptotic behaviors 
given in the text (|48j |49]) . 
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